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ABSTRACT 

We assess the impact of starspots on the evolution of late-type stars during the pre-main sequence 
(pre-MS) using a modified stellar evolution code. We find that heavily spotted models of mass 0.1- 
i.2M0 are inflated by up to 10% during the pre-MS, and up to 4% and 9% for fully- and partially- 
convective stars at the zero-age MS, consistent with measurements from active eclipsing binary sys¬ 
tems. Spots similarly decrease stellar luminosity and Teg, causing isochrone-derived masses to be 
under-estimated by up to a factor of 2x, and ages to be under-estimated by a factor of 2-lOx, at 
3 Myr. Consequently, pre-MS clusters and their active stars are systematically older and more massive 
than often reported. Cluster ages derived with the lithium depletion boundary technique are erro¬ 
neously young by ^ 15% and 10% at 30 and 100 Myr respectively, if 50% spotted stars are interpreted 
with un-spotted models. Finally, lithium depletion is suppressed in spotted stars with radiative cores, 
leading to a fixed-temperature lithium dispersion on the MS if a range of spot properties are present 
on the pre-MS. Such dispersions are large enough to explain Li abundance spreads seen in young 
open clusters, and imply a range of radii at fixed mass and age during the pre-MS Li burning epoch. 

By extension, this implies that mass, composition, and age do not uniquely specify the HR diagram 
location of pre-MS stars. 


1. INTRODUCTION 

Precision measurements of fundamental stellar param¬ 
eters active, late-type stars are systematically larger by 
~5-15% than predictions from modern evolution codes, 
a phenomenon we refer to as radius inflation. This ef- 


feet has been claimed in open clusters (e.g. 

Jackson & 

Jeffries||2014aD and in field stars (e.g. Popper 

19971 'i'or- 

res & Ribas||2UL)21 |Ribas||2003l ILopez-Morafes & Ribas 

2005 Morales et al. 2008| Boyajian et al. 

20121): on the 

main sequence (MS), and the pre-MS (e.g. 

Stassun et al. 


20071; and for s tars above and below the fulfy-convective 


bound ary (e.g. Torres et al. 2010 Feiden & Chaboyer 


2012a I. The cause remains controversial, but a vital clue 


comes from the claimed correlation of radius inflation 


with magnetic activity (Torres et al. 


2007 Morales et al. 


Lopez -Morales 
20081 IStassun et al.||2012|). Mech¬ 


anisms proposed to induce this correlat ion include the 


dence is accruing that such radius anomalies are real, in¬ 
terpreters must be caut ious as complication s such as mis¬ 


estimated opacity (e.g. Terrien et al. , unaccounted 


systematics in interferometic measuremen 


s (Casagrande 


et al.||2014|, and the impact of spots on radius m easure- 


mist aken for signatures of radius in flation. 


In Somers & Pinsonneault (2014, SP14 hereafter), we 
found that radius inflation during the pre-MS has a 
strong impact lithium abundances of young stars. Ra¬ 
dius inflation reduces the temperature at the base of the 
surface convection zone during the pre-MS, leading to 
a suppressed rate of Li destruction in standard stellar 


models. If stars of equal mass and age have different in¬ 
trinsic sizes during the pre-MS, they will destroy differ¬ 
ent quantities of Li during the pre-MS, thus providing a 
mechanism for inducing the Li abundance dispersion ob- 


clusters (e.g. 

Soderblom et al.||1993 

sonneault 

(2015 

SP15 hereafter), we 


tation al mixing on the pre-MS (e.g. Pinsonneault et al. 


19901. This implies that if rotation correlates with an in¬ 


lardo & Baraffe 

2007 

), magnetic fields adding internal 

pressure support (e.g. IMulfan &: MacDonald 

2001 

), and 

starspots (e.g 

Torres et af.||2006p, though eac 

n are 

likely 


creas ed radius on the pre-MS, as claimed in som e works 
(e.g. Littlefair et al. 2011 Cottaar et al. 2014), the Li- 
rotation correlation observed in the Pleiades could also 
be produced. This is important because rotational mix¬ 
ing is correlated with rotation rate, producing a trend of 
opposite sign relative to the data. Both of our previous 
papers were constructed as general investigations of the 
effect of radius inflation on Li abundances, conducted 


by varying the mixing len gth parameter (e.g. Chabrier 


ments m eclipsing binaries ((Morales et al. 112010) may be ( Gough &: Tayler 

. r. . r. 1—1^-1—II— -II-1 hppn arcrnArl fn at 


Gallardo & Baraffe 2007), and as such did not model a 
unique physical mechanism. In this paper, we expand 
upon this work by exploring a specific potential cause of 
radius inflation: starspots. 

Spots are the visual manifestation of concentrated 
magnetic fields in the convective envelopes of late-type 
stars. These fields inhibit the flow of energy in specific 
regions near the surface, inducing a st ructural response 


radii of stars (e.g. Spruit 

1982a 

Spruit & Weiss 

1986 

I, 

but in only a few cases have their effect on the pre-MS 

been examined (e.g. Jackson & Jeffries 

2014a 

I- 

Spruit 


& Weiss (1986) proposed a method for introducing the 


physical effects of spots into stellar evolution codes, with 
the goal of exploring the colors of stars with black spots 
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at the zero-age MS (ZAMS). In this paper, we imple¬ 
ment a similar treatment into our evolutionary code, and 
explore the impact of spottiness on pre-MS and ZAMS 
stellar parameters, the inference of masses and ages from 
the Hertzsprung-Russell (HR) diagram, the observed col¬ 
ors of stars with non-zero spot temperatures, and the 
destruction of lithium. By comparing our model predic¬ 
tions to data, a consistent picture emerges: pre-MS stars 
have a range of radii at fixed mass and age, correlated 
with rotation, which persists at the ZAMS. 

The text is organized as follows. We discuss our evo¬ 
lutionary models in with a particular emphasis on 
our treatment of spots ( p.2[ ). With our models calcu¬ 
lated, we hrst consider the impact of spots during the 
pre-MS and at the ZAMS on stellar radiiMuminosities, 
and TeffS in the mass range O.I-I. 2 M 0 (®. Next, we 
discuss how spots might introduce errors into the q uan- 
tification of masses and ages using isoc hrones ( |4.1[ ), and 
how spots affect observed colors ( ^4.2[ ). Then, we exam¬ 
ine how spots impact the destruction of lithium, both 
in the context of age determinations with the lithium 
depletion boundary technique, and in the context of Li 
abundance spreads observed in young MS clusters ((J^. 
In we synthesize the results of (|3]--[^ compare our 
quantitative findings to other work in th^iterature, and 
discuss some outstanding questions for starspot model¬ 
ing. Finally, we summarize our conclusions in © 


star will begin to contract towards the hydrogen burning 
main sequence. Our zero point provides a consistent cri¬ 
teria for the starting age of all models, and ensures that 
models of equal mass but differing spot properties are in 
equivalent evolutionary states when compared at early 
ages. 


2.2. Spots 

2.2.1. Implementation 

Our spot me t hodol ogy, based on the treatment by 
Spruit & Weiss (1986), incorporates two effects: the in- 
homogeneous photosphere in the presence of spots, and 
the inhibition of convection in magnetic regions. To sim¬ 
ulate this behavior, we divide our stellar models into two 
zones: un-spotted and spotted. We account for the flux 
redistribution by explicitly altering the radiative trans¬ 
port gradient in the surface convection zones of our mod¬ 
els, and the inhomogeneous surface by properly treating 
the outer boundary condition (see below). 

The spots in our models are controlled by two param¬ 
eters: the filling factor fspot, defined as the ratio of the 
spotted surface to the total surface area, and the spot 
temperature contrast ratio Xgpotj defined as the tempera¬ 
ture of the spotted surface divided by the temperature of 
the un-spotted photosphere. The surface-weighted aver¬ 
age of the two zone temperatures defines the temperature 
of the model at a given radius. 


2. METHODS 

2.1. Standard stellar models 

We calculate non-rotating, solar-metallicity stellar 
mod els with the Yale Rotating Evolutio n Code (YREC; 


e.g. van Saders & Pinsonneault 2013). For the pur¬ 


poses or"caIibratIbn7”we~assumetne~Sun is unaffected 
by the presence of spots; this is reasonable, given the 
Sun’s minimal spot coverage. We adopt a present day 
solar heavy element abun dance of ZjX = 0.02292 from 
Grevesse & Sauval (19981, and choose the hydrogen mass 
traction (A ) and the mixing length coefficient (oml) that 
reproduce the solar luminosity and radius for a IMq 
model at 4.568 Gyr. The final solar calibrated values 
are A = 0.71172, Z = 0.018923, and a© = 1.925 31. Our 
models use the 2006 OPAL equation of state (Rogers, 


1996 


Rogers fc Nayfonov 


Tom|AlTard et all 1 1997 T hign tem- 


imr, 


Swenson & Iglesias 
model atmospheres 
perature opacities from the opacity proje ct (Mendoza et 
al.|200 7|) , low temperature opacities fro m|h'erguson et al. 
(20051), and nuclear cross-sections from Adelberger et al. 


. We tre at electron screening using the method of 
Saipeterj (|1954|). 

Uur models are initialized above the Hayashi track, 
and evolved to the age at which 99% of t he initial deu - 
terium abundance has been destroyed (e.g.|Stahler|1988|). 
We define this time as t = 0 for each individual model. 
This approach represents a physical scenario: as proto¬ 
stars contract and build up mass through accretion, their 
interiors increase in temperature. Eventually, the core 
becomes hot enough to fuse deuterium, a process which 
partially supports the star against gravitational contrac¬ 
tion; this is the so-called “deuterium birth line”. Proto¬ 
stars travel along this sequence towards higher mass as 
accretion continues to supply fresh deuterium fuel to the 
central engine. Once the accretion rate subsides, the 


T" = (1 - fspot)TLb + UotT^pou (1) 

Tspot is the area-averaged temperature of the prenum- 
bral and umbral spot regions, and Tamb is the tempera¬ 
ture of the ambient, un-spotted region. By substituting 
Tspot = XspotTamb, we find an expression for the total 
luminosity. 


L — AttR asBTg^^ij{l fspot + fspotXgpgf), (2) 
and the flux in the ambient regions. 


Tamb — TiQi/il fspot T fspotXspot)' 


(3) 


This formulation implies a modeling degeneracy be¬ 
tween fspot and Xspot, since the flux redistribution is ac¬ 
tually controlled by the redistribution parameter, 


C^spot — (1 fspot + fspotXgpot)- (4) 

otspot can be thought of as one minus the effective 
blocking area of totally black {Xspot = 0) spots, a pa¬ 
rameter typically called /3 (= 1 — Uspot)- Eq. [^governs 
the flux enhancement in the non-magnetic regions of the 
model. We implement this enhancement into our evolu¬ 
tion code by dividing the radiative gradient in convective 
regions by Ospot, 


^ rad,amb — ^ rad/Otspot- 


(5) 


Since fspot and Xspot fall in the range O-l by defini¬ 
tion, aspot is always less than or equal to one, meaning 
that spots can only increase the radiative flux in the un¬ 
spotted regions. The luminosity and effective tempera¬ 
ture of the star can be found by solving Eqs. and 
using the values of Tamb, Tspot, and Xspot at the surface 
where the optical depth t = 2/2,. To model the changing 
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temperature of spotted regions with depth, we hold the 
difference in temperature between the two zones fixed 
(Tamb — Tspot = const.). The result is that spots are only 
signihcant near the surface, where T is low. We tes t the 
importance of spot depth on our models in §6.3.1[ and 
find that assumptions about the subsurface morphology 
impact radius, luminosity, and Teff predictions by ^ 1 % 
at most. 

Our spot treatment involves a crucial implicit assump¬ 
tion: the spotted and un-spotted atmospheric regions 
are in pressure equilibrium. This means that if the spot¬ 
ted regions contain a magnetic field which produces a 
pressure Pspot,mag, the gas pressure in the two regions 
will be such that Pamb.gas — Pspot,gas T .Pspo^,maJ~l This 
requirement permits us to obtain the surface pressure in 
the un-spotted region, and apply it equally in both zones. 
Psurf is obtained in our evolution code by using a model 
atmosphere, which returns a pressure Psurf for a given 
Teff and logg. For consistency, we call our atmosphere 
routine using the temperature of the un-spotted, ambient 
surface: 


Psurf — Psurf (Pamb,surf 1 fogff)- 


( 6 ) 


Our approach differs from previous studies which have 
considered the effects of magneti c activity on stellar 
properties. | Spruit & Weiss | (|1986|) used a similar spot 
formulation to investigate the impact on stellar param¬ 
eters, but focused primarily on the ZAMS, and did not 
consider the effect of non-zero spot temperatures on the 
location of stars in the color-magnitude diagram. In a 
recent series of papers, |Feiden & Chaboyer| (2012b|) ex¬ 
amined the consequences ot magnetic helds in stellar in¬ 
teriors by directly incorporating a 1-dimensional approxi- 
mation of the induction equation into their model physics 
(Feiden & Chaboyer 2013 20141. This approach differs 
from ours by considering a one-zone, homogeneous sur¬ 
face boundary condition, compared to the two-zone, in¬ 
homogeneous approach we have adopted. However, since 
magnetic fields are the underlying cause of surface spots, 
our work compliments theirs by approaching the same 
mechanism from a differen t perspective. In addition, 
Jackson & Jeffries (2014a|b) have recently examined the 
impact ot surface spots on the pre-MS radii of late-type 
stars, and the lithium destruction rates of stars with M 
< O.5M0. Their procedure relied on calibrating poly¬ 
tropic models on published evolutionary tracks, and then 
considering the impact of flux-blocking spots at the sur¬ 
face. By contrast, we employ a self consistent evolution¬ 
ary code to produce our calc ulati ons, and thus anticipate 
somewhat different results (j |6.2[ ). Each of these studies 
differed in their treatment from our present study, and 
we consider our results complimentary to theirs. 


2.2.2. Adopted spot properties 

In this section, we select a range of spot properties to 
investigate in our models. We are interested in the range 
of stellar properties that could be produced by spots, 
as well as their dependence on evolutionary state and 
mass. Most measured filling factors to date have been 

^ A further assumption of magnetic and gas equiparition there¬ 
fore implies Pamb gas = 2 X Pspot gas, a fact we make use of in 


for dwarfs , sub-giants, and giants (see tab. 9 in Berdyug-] 


ma 


20051. While very few pre-MS stars have similar 


measurements, perhaps the most direct comparison are 
the spotted, tidally-synchronized sub-giants known as RS 
CVn stars. RS CVns are structurally similar to pre-MS 
stars in the outer regions, as they have deep convective 
envelopes and low gravities relative to the MS, and thus 
host similar environments for spot emergence. Claimed 
filling factor among this class reach 40 to 5 0%, such as in 


the cases of HU Vir, UX Ari, and HR1099 (O’Neal et al. 


1998 2001), demonstrating that deep envOopes paired 


with moderate (5 — 10 d) rotation periods can produce a 
high level of spot coverage. The T Tauri star V4 10 Tan 


has a measu red surface spot coverage of 32-41% (Petrov 


et al. 1994), supporting the notion that pre-MS stars 


can be similarly spotted. Furthermore, Zeeman-Doppler 
imaging has confirmed avera ge surface hel ds of up to a 
few kG on pre-MS stars (e.g. Reinersj 20121, comparable 
to the surface fields of highly active, heavily spotted MS 
stars (see below). Collectively, these arguments support 
the hypothesis that pre-MS stars can host spots which 
approach 50% coverage. 

Similar filling factors have been claimed for young MS 
stars. For example, the rapidly rotating (Trot 1 to 
4 d) dwarfs LQ Hya, V83 3 Tan, and EQ Ve r host fill¬ 
ing factors of 40 to 50 % (jO’Neal et al.||2001| |2004|). It 
may appear surprising at hrst that the slower rotating 
RS CVn sta rs (e.g. Trot ~ 10.4 d for HU Vir; |0 Neal 


et al. 1998) can be as spotted as the faster rotating 
ZAlVlS stars. However, the convective overturn timescale 
of sub-giants is much longer than MS dwarfs, implying a 
smaller Rossby number (ratio of rotation rate to overturn 
timescale) at fixed rotation period. Magnetic properties 
appear to scale more closel y with Rossby number than 
with absolu te rotation (e.g. Noyes et al.||1984' Pizzolato 


er£l][2003t , so it is reasonable to expect equivalent spot 


properties at slower rotation rates for sub-giants. This 
argument applies equally well to pre-MS stars, whose 
average rot ation rates are slower at 3 Myr than at the 
ZAMS (e.g. Littlefair et al.pOlO ), but whose convection 
zones are deeper. 

We conclude that fspot = 0.5 is a reasonable assump¬ 
tion for the upper limit of pre-MS and ZAMS filling fac¬ 
tors, and we adopt this. We further adopt Xspot = 0.8, 
characteristic of the range of observed spot te mperature 
ratios (0.65 to 0.85; e.g. |Berdyugina| |2005|) . A fuller 
treatment would include a mapping of spot tempera¬ 
ture and filling factor as a function of age, metallicity, 
mass, and rotation rate; the current exercise is aimed 
at a different purpose, namely constraining the potential 
impact of star spots on observables. We do note that 
the Pleiades lithium dispersion (see does provide ev¬ 
idence for a range o f spot properties m pre-MS stars, as 
we argued in|SP15[ 


3. EFFECT ON STELLAR RADII 
3.1. The pre-main sequence 

We first explore the impact of spots on stellar model 
radii during the pre-MS. Evolutionary tracks represent¬ 
ing both fully convective stars (M < O.dSM©) and stars 
which develop radiative cores (O.35M0< M < I. 2 M 0 ) 
are shown in the upper panel of frame A in Fig. Solid 
black lines represent the radius as a function of time for 
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log(Age/Myr) 




Figure 1. The impact of spots on radii, luminosities, and Tgfjs of 
stars in the mass range O.1-1.2M0. A) The time evolution of the 
radius of un-spotted (solid) and spotted (dashed) models of various 
masses are shown in the top. The bottom shows evolution of the 
difference in radius between the spotted and un-spotted models for 
fully convective (dotted) and partially convective (solid) models. 
AR steadily grows during the early pre-MS, before stabilizing at 
the entrance to the MS. B) The top shows zero-age MS radius 
anomalies across our mass range, for three different values of /spot- 
The bottom plots the luminosity and anomalies resulting from 
the case of fspot = 0.5, respectively shown in blue and red. C) The 
mass-radius relationship for un-spotted (light) and spotted (dark) 
models. The bottom panel shows that the mass-dependence of the 
structure of the relationship is affected by spots; this leads to the 
strong mass dependence of AR shown in panel B. 


un-spotted models, and the dashed black lines represent 
the same for models half covered in 80% temperature 
contrast spots i^fspot — 0.5, Xgpot — 0.8, o^spot ^ O.T), 
referred to hereafter as our base case. The difference in 
radius between equal mass, un-spotted and spotted mod¬ 
els as a function of age is reported in the lower panel. 

Spots lead to an inflated radius at all ages, but the 
magnitude evolves during the pre-MS. For each mass, 
the radius enhancement gradually increases with age, ini¬ 
tially peaking around 8 - 12 %, and finally stabilizing on 
the MS. The age dependence of pre-MS radius inflation 
results primarily from a timing argument: the rate of pre- 
MS contraction is proportional to stellar luminosity, and 
because un-spotted stars are more luminous than their 
spotted counterparts on the pre-MS (see below), they will 
have contracted further at a given age. This explains the 
gradual increase in AR with time, and leads to a conspic¬ 
uous transient spike in AR at the onset of MS burning, 
as the un-spotted star has reached the ZAMS while the 
spotted star remains on the late pre-MS. However the in¬ 
ternal structure is also altered in a mass dependent fash¬ 
ion: changes to t he efficiency of convection impact higher 
mass models (e.g. Andronov & Pinsonneault 2004), while 
changes to the surf ace boundary conditio n are important 
at lower mass (e.g. [Spruit fc Weiss 1986). 

Spotted models show corresponding anomalies in both 
Teff and luminosity. The average Teg for our base case is 
lower by up to 7% for all masses, and corresponds to a 
zero point shift in the Hayashi track. Pre-MS luminosity 
anomalies therefore evolve according to L/R^ ^ const. 
Luminosities of the spotted models are lower by up to 
18% and 10% during the early and late pre-MS for high 
mass stars, and lower by up to 25% and 10% during the 
early and late pre-MS for low mass stars. We explore 
consequences of these Teff and L anomalies in Q and ^ 

3.2. The zero-age main sequence 

Once at the ZAMS, radius anomalies converge to a sta¬ 
ble, mass-dependent pattern. The black line in the upper 
panel of frame B in Fig. shows our base case radius 
anomalies at the age when 0 . 1 % of the core hydrogen has 
been burned. As can be seen, fully convective stars have 
a weakly mass-dependent anomaly of ^4%, and higher 
mass stars generally have larger inflation factors with lo¬ 
cal maxima at O. 5 M 0 and I. 2 M 0 , and a local minimum 
at O. 7 M 0 . Anomalies resulting from different values of 
fspot are also shown. Unsurprisingly, the radius enhance¬ 
ment increases self-similarly with greater spot coverage, 
and the change in Ai? with fspot is approximately lin- 
ear. These r esults are generally consistent with [Spruit fc| 
Weiss (1986), though the second minimum around O. 8 M 0 
is not present in their ZAMS models. This may be due 
to differences in the adopted standard model physics. 

This radius enhancement in spotted models is predom¬ 
inately a result of the altered boundary conditions at the 
surface. The photospheric pressure is higher for the hot¬ 
ter ambient photosphere than it would be for a model 
with Tsurf = Teff at fixed radius. The model expands 
as a result in order to reach equilibrium (P oc M'^/R'^ 
in homology), which lies at higher radius and lower sur¬ 
face pressure. This reduced P leads to a corresponding 
decrease in P throughout the surface convective zone, 
due to its adiabatic stratification. For fully convective 
stars, the pressure reduction reaches the nuclear burn- 
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Figure 2. Left: The mass-radius relationship for un-spotted (black) and spotted (blue) models at the ZAMS, and the 10 Gyr relation for 
un-spotted models (red dashed). These are compared to data from eclipsing binaries (see the text for references). Right: The same as the 
left, but detrended against the ZAMS, un-spotted isochrone. Dotted purple represents the sum of age and spot effects. Two-thirds of the 
data lie within the spotted models, and all but one lie within the combined effects of activity and age. 


ing core where it directly suppresses the rate of reac¬ 
tions, and thus lowers the luminosity. The blue line in 
the lower panel of frame B shows AL as a function of 
mass, revealing a significant decrease in the luminosity 
of fully convective stars due to our base case spots. This 
effect is less pronounced in higher mass stars, whose con¬ 
vection zones have receded from the core and therefore 


exert less influen ce on the central engine (Kippenhahn 
& Weigert |1990 ). The decreasing luminosity anomaly 


from U.d-U.bJVlQ leads to a spike in the radius and tem¬ 
perature (red line in frame B) anomalies, which peak 
around O.SMq and O.TMq respectively. Note that above 
^O. 5 M 0 , AR « —2 Ar, in good ayeem ent with obser¬ 
vational studies (e.g. |Morales et al.|2008|). 

The peculiar structure of the dependence of AR on 
mass can be understood by considering the spotted and 
un-spotted mass-radius relationships plotted in the up¬ 
per panel of frame C in Fig. As can be seen, between 


0.1 and I.OM 0 the slope of the mass-radius relation at 
first decreases, then increases, decreases, and increases 
again. These “wiggles” result from mass-dependent 
changes in stellar structure and energy generation, such 
as the onset of a radiative core, the decreasing mass of the 
convective envelope, and the burgeoning dominance of 
CNO hydrogen processing. The lower panel shows both 
mass-radius relations detrended with the line M = i? (in 
solar units). The inclusion of spots alters the masses at 
which these transitions occur, shifting the onset of slope 
changes to lower mass, and thus leading to significant 
structure in the mass dependence of radius inflation. 


3.3. Comparison with data 

In this section, we compare the mass-radius relation¬ 
ship predicted by our models to literature data. The 
left panel of Fig. displays our predicted relationships 
for un-spotted (black line) and spotted (blue line) mod¬ 
els at the ZAMS, and for un-spotted models at 10 Gyr 


(dashed red). These ar e compared to M < l .OMfn eclips- 
ing binary data from iTorres et al.l (|2010|) and iFeiden 


fc Chaboyerj (2012a) (red and black, respectively). The 


binaries generally form a tight sequence bracketing the 
spotted isochrone below O. 8 M 0 , and fan out into a larger 
range above O. 8 M 0 , attesting the importance of age ef¬ 
fects at a solar mass and above. In the right panel, 
the data falling within the dotted parallelogram are de¬ 
trended with the standard ZAMS isochrone, and com¬ 
pared to the spotted and 10 Gyr isochrones. The data 
occupy an envelope ^ 4 — 7% at the low mass end, and 
increase to ~ 12 % at O. 7 M 0 , demonstrating the inflated 
radius problem discussed in Two-thirds of the de¬ 
trended stars fall within the shaded blue region, and all 
but one fall within the combined effects of age and spot^ 
(dotted purple line). This shows that the radius anoma¬ 
lies resulting from realistic spot parameters in our models 
are consistent with empirical constraints on radius infla¬ 
tion. 

We find no clear evidence for higher av erag e anomalies 
around O. 5 M 0 compared to O. 75 M 0 (see |3.2| , but a host 
of systematic issues forbid any present conclusions. For 
example, stars of equal mass but unequal metal abun¬ 
dance will have differing radii, but we And that this dis¬ 
persion is only ~ 2% for [Fe/H] S [—0.3, -1-0.3]. Addition¬ 
ally, surface spots can increase the inferred transit time 
for eclipsing binaries, thus leading to syste matic over¬ 


est imate s of up to ^3% in the radius (e.g. Morales et 


al. 2010 1 . Finally, age effects will contribute to inferred 


anomal ies, as suggested by Fig. [Feiden ik Chabo}^ 
( 2012 a) have recently argued that accounting for these 
effects reduces most radius anomalies to around 5%, in¬ 
stead of the often quoted range of 5-15%, in good agree¬ 
ment with our expectations. In any case, it appears that 


^ Old stars are normally not expected to be heavily spotted, but 
short period eclipsing binaries are often tidally locked, permitting 
rapid rotation at late ages. 
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Figure 3. The impact of spots on the HR diagram. Left: Un-spotted (red) and spotted (blue) evolutionary tracks and 3 Myr isochrones 
for a range of masses. Spots displace stars to the lower right (cooler and less luminous) at fixed age. As a result, a dispersion in luminosity 
at fixed Teff develops across the mass function, and spotted stars appear less massive and younger when interpreted with un-spotted 
isochrones. Right: Correction factors which describe the errors incurred when measuring the mass and age of spotted stars with un-spotted 
isochrones. At a given mass (x-axis) and age (red/orange/yellow) inferred from un-spotted isochrones, the y-axis reports the necessary 
corrections. For example, a heavily spotted star which appears to have a mass of O.SMq and an age of 3 Myr will have respective corrections 
(green crosses) of 1.5 (implying a true mass of 1.5 x 0.5 Mq ~ 0.75 Mq), and 2.7 (implying a true age of 2.7 x 3 Myr ~ 8 Myr). 


spot-induced radius effects are of a similar magnitude 
as the well-known radius anomalies in eclipsing binaries, 
and should not be neglected when interpreting these ob¬ 
servations. 


4. EFFECT ON INFERRED MASSES, AGES, AND COLORS 

Measuring the masses and ages of pre-MS stars is 
a challenging endeavor, given the rapid evolut ionary 
timescales a nd a host of observational issues (e.g. Hart- 
mann]2doT ), but is of great importance for studies of the 
initial mass function, circum-stellar discs, the early an¬ 
gular momentum evolution of stars, and the formation of 
planets. The ideal approach would be to compare cali¬ 
brated pre-MS isochrones with the HR diagram of a clus¬ 
ter, and to read off the masses and ages of its members. 
In reality, this procedure is complicated by a variety of 
observational and theoretical issues, which broadly fall 
into three categories. 

Observational errors. In young clusters, extinction is 
highly variable from source to source, binarity status is 
often unknown, and stellar magnitudes may vary on short 
and long timescales due to rotation, disc obscuration, and 
accretion. These effects form a base level of uncertainty. 

Theoretical issues. Pre-MS evolutionary tracks lack a 
good calibration data set, and different suites make vary¬ 
ing predictions for the birth line location, the bolometric 
corrections, and the rate of Hayashi contraction. This 
uncertainty severely limits one’s ability to derive masses 
and ages in young clusters. 

Non-uniform parameters. Even if these issues were re¬ 
solved, the possibility remains that stellar parameters 
such as luminosity and radius vary between cluster mem¬ 
bers at fixed mass and age, perhaps due to physics asso¬ 
ciated with rotation, magnetic activity, or accretion. 

The picture is further complicated by intra-cluster age 


spreads, which mimic the appearance of non-uniform 
stellar parameters by producing luminosity dispersions 
at fixed Teg. In order to separate these effects, it is im¬ 
portant to asses the impact that non-uniform parameters 
would have on coeval stars. 

In ^ we showed that spots induce shifts to the radii, 
luminosities and T^ss of stars during the pre-MS, and can 
therefore generate non-uniform parameters among young 
cluster members. In this section, we explore the im¬ 
pact of these shifts on the Hertzprung-Russell and color- 
magnitude diagrams of young clusters, and discuss their 
implications for understanding pre-MS stars. 


4.1. Spot effects in the Hertzsprung-Russell diagram 

We first constrain ourselves to the theoretical HR dia¬ 
gram, and consider the impact of spots on the tempera¬ 
ture and luminosity of our stellar models. This for now 
ignores the impediments to obt aining these parameters 
brought on by spots (see (4.2), and assumes the true 


surface and L are known. 

The red dashed lines in the left panel of Fig. show 
un-spotted tracks for a range of masses, run from the deu¬ 
terium birth line (dotted red line) to 50 Myr. The blue 
dashed lines represent the same for our base case spot¬ 
ted models. Isochrones of 3 Myr are over-plotted in solid 
red and blue, with masses periodically enumerated. The 
results of (j^are reflected in this figure: spots reduce the 
intrinsic luminosity and Teff of stars in a mass-dependent 
fashion. The fractional displacement ranges from 
^6-8% from the highest to the lowest mass, which corre¬ 
sponds to a nearly mass-independent decrement of 270 K. 
The fractional luminosity reduction ranges from ~10- 
22% from high to low mass, corresponding to shifts of 
0.1 - O.OI Lq. A range of spot filling factors will there¬ 
fore induce a significant dispersion in L at fixed T^g. 
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However, even in our extreme base case the dispersion 
does not reach the order-of- magnitude spread o bserved 
in some young clusters (e.g. Hillenbrand 19971. Thus 
spots may be an important contributor to young cluster 
dispersions, but not their sole cause. We emphasize that 
these shifts are not an observational bias, but a change 
in the true temperature and luminosity of pre-MS stars. 

Pre-MS Hayashi evolution proceeds at approximately 
fixed Teg, with L declining with age. Our models demon¬ 
strate this behavior, but also reveal an important result: 
spots alter the Tgg at which this evolution occurs. Spots 
displace stars downward and to the right in the plane 
of the HR diagram, and consequently spotted stars ap¬ 
pear younger and less massive when interpreted with 
standard isochrones. To quantify these errors, we use 
the un-spotted evolutionary tracks of Fig. [^to infer the 
masses and ages of a range base case spotted^models, and 
derive a set of correction factors to the resulting param¬ 
eters. The results appear in the right panel of Fig. 
These factors operate in the following sense: if the mass 
and age of a 50% spotted star are derived from its HR 
diagram location using un-spotted tracks, the true mass 
and age will be given by the derived values multiplied by 
their respective correction factors (see the caption for a 
specific example). Factors for a variety of masses, ages, 
and filling factors can be found in Table 

Mass correction factors (solid lines) for solar-type stars 
are as large as 1.3x, but decrease with inferred age 
since spotted and un-spotted evolutionary tracks begin 
to overlap near the horizontal Henyey track (e.g. the 
O.9M0 tracks). Mass errors increase substantially to¬ 
wards the fully convective regime, where spotted star 
masses can be under-estimated by up to 2.2 x, and never 
less than 1.7x, regardless of age. The larger errors for 
low mass stars result from the stronger impact of spots on 
fully convective luminosities, and the temporal stability 
of the errors results from the long contraction timescale 
of low mass stars. Age correction factors (dashed lines) 
are more severe. For solar-type stars, ages are under¬ 
estimated by up to 2.5X, and become more accurate as 
stars get older. At O.2M0, ages will be under-estimated 
by up to a factor of 10, with the largest errors coming at 
the earliest epochs. This order-of-magnitude shift results 
from the large fractional errors incurred when consider¬ 
ing values near zero - for instance, an under-estimate 
of 1 Myr for a 1.1 Myr old star produces a correction 
factor of 11. They also relate to the decreasing luminos¬ 
ity of the deuterium birth line with mass: spotted stars 
are born at a higher luminosity at hxed Tea, and thus 
lie above un-spotted stars at fixed age (Fig. [^. Conse¬ 
quently, spotted stars interpreted with un-spotted tracks 
will appear far younger than their true ages, thus requir¬ 
ing a correction factor much larger than unity. We note 
that the magnitude of the age corrections depends some¬ 
what on how one defines t = 0, but that mass correction 
factors are largely insensitive to this choice. 

We note an interesting prediction of these models. In 
our physical scenario, proto-stars travel up the deuterium 
birth-line during the assembly phase, slowly gaining mass 
through accretion. They are prevented from substan¬ 
tial contraction by deuterium burning feedback, and 
prevented from persisting in a higher luminos ity state 


by the short Kelvin-Helmh oltz time-scales (e.g. Baraffe, 


Chabrier & Gallardo 2009). The region of the HR dia- 


Table 1 

Mass and Age Correction Factors 


Inferred Mass 
(Mr) 

Mass Correction 
Factors 

Age Correction 
Factors 

fspot — 1/6 

1 Myr 3 Myr 10 Myr 

1 Myr 3 Myr 10 Myr 

0.20 

1.263 

1.275 

1.296 

3.589 

2.311 

1.881 

0.22 

1.249 

1.263 

1.289 

3.280 

2.191 

1.835 

0.24 

1.234 

1.251 

1.281 

3.010 

2.084 

1.790 

0.26 

1.220 

1.239 

1.270 

2.777 

1.986 

1.740 

0.28 

1.205 

1.226 

1.263 

2.561 

1.893 

1.704 

0.30 

1.194 

1.214 

1.255 

2.397 

1.818 

1.668 

0.35 

1.175 

1.194 

1.236 

2.120 

1.686 

1.597 

0.40 

1.163 

1.181 

1.217 

1.939 

1.600 

1.535 

0.45 

1.151 

1.172 

1.198 

1.782 

1.532 

1.495 

0.50 

1.137 

1.160 

1.178 

1.648 

1.464 

1.449 

0.55 

1.123 

1.148 

1.155 

1.535 

1.407 

1.406 

0.60 

1.106 

1.135 

1.133 

1.434 

1.357 

1.368 

0.65 

1.095 

1.119 

1.111 

1.383 

1.312 

1.332 

0.70 

1.093 

1.109 

1.090 

1.360 

1.286 

1.295 

0.80 

1.093 

1.097 

1.059 

1.335 

1.264 

1.254 

0.90 

1.092 

1.085 

1.036 

1.313 

1.242 

1.196 

1.00 

1.089 

1.074 

1.013 

1.294 

1.226 

1.163 

1.10 

1.088 

1.063 

1.002 

1.286 

1.208 

1.135 

1.20 


1.052 

0.994 


1.193 

1.127 

fspot — 1/6 

1 Myr 3 Myr 10 Myr 

1 Myr 3 Myr 10 Myr 

0.20 

1.593 

1.633 

1.710 

7.349 

4.321 

3.419 

0.22 

1.548 

1.595 

1.681 

6.412 

3.937 

3.244 

0.24 

1.510 

1.560 

1.650 

5.686 

3.622 

3.083 

0.26 

1.477 

1.529 

1.620 

5.104 

3.358 

2.938 

0.28 

1.450 

1.501 

1.592 

4.641 

3.139 

2.813 

0.30 

1.427 

1.477 

1.564 

4.268 

2.958 

2.702 

0.35 

1.388 

1.435 

1.499 

3.606 

2.635 

2.487 

0.40 

1.355 

1.401 

1.434 

3.114 

2.399 

2.291 

0.45 

1.321 

1.370 

1.375 

2.714 

2.204 

2.132 

0.50 

1.284 

1.338 

1.321 

2.383 

2.034 

1.999 

0.55 

1.248 

1.303 

1.271 

2.125 

1.887 

1.883 

0.60 

1.220 

1.271 

1.224 

1.949 

1.772 

1.781 

0.65 

1.206 

1.243 

1.184 

1.856 

1.689 

1.694 

0.70 

1.202 

1.223 

1.148 

1.801 

1.639 

1.627 

0.80 

1.195 

1.191 

1.092 

1.723 

1.571 

1.524 

0.90 

1.183 

1.161 

1.046 

1.660 

1.513 

1.403 

1.00 

1.182 

1.132 

1.017 

1.644 

1.460 

1.327 

1.10 

1.182 

1.110 

1.000 

1.640 

1.428 

1.280 

1.20 



0.984 



1.279 

fspot = 1/2 

1 Myr 3 Myr 10 Myr 

1 Myr 3 Myr 10 Myr 

0.20 

2.005 

2.093 

2.156 

12.819 

7.453 

5.968 

0.22 

1.931 

2.023 

2.074 

11.039 

6.649 

5.490 

0.24 

1.866 

1.961 

2.000 

9.639 

5.993 

5.076 

0.26 

1.812 

1.904 

1.931 

8.526 

5.453 

4.721 

0.28 

1.765 

1.854 

1.867 

7.630 

5.003 

4.413 

0.30 

1.725 

1.809 

1.808 

6.888 

4.627 

4.139 

0.35 

1.643 

1.719 

1.680 

5.523 

3.928 

3.591 

0.40 

1.573 

1.644 

1.571 

4.545 

3.419 

3.177 

0.45 

1.502 

1.575 

1.478 

3.789 

3.017 

2.864 

0.50 

1.439 

1.510 

1.401 

3.229 

2.692 

2.609 

0.55 

1.390 

1.454 

1.333 

2.841 

2.450 

2.407 

0.60 

1.352 

1.405 

1.275 

2.571 

2.270 

2.245 

0.65 

1.329 

1.362 

1.224 

2.409 

2.134 

2.110 

0.70 

1.313 

1.325 

1.179 

2.291 

2.038 

1.992 

0.80 

1.286 

1.265 

1.105 

2.134 

1.893 

1.790 

0.90 

1.285 

1.211 

1.051 

2.108 

1.772 

1.612 

1.00 

1.286 

1.177 

1.018 

2.108 

1.725 

1.491 

1.10 


1.150 

0.994 


1.703 

1.454 

1.20 



0.968 



1.480 


Note. — Mass and age correction factors for three spot fill¬ 
ing factors at vari ous inferred masses and ages. Correct useage 
is detailed in ^4.1| Ellipses indicate inferred mass and age points 
outside of the interpolation grid. 






















Figure 4. The impact of spots on the color-magnitude diagram. Using empirical bolometric corrections, our un-spotted (red) models are 
directly transformed into three different color coordinates using a two-zone spot model. In color coordinates, the spotted models do not 
appear as cool due to the contribution to short wavelength emission from the un-spotted photosphere. The spot models appear bluer when 
shorter wave-lengths are used, due to the progressively stronger sampling of the Wein tail of the blackbody function. The bottom right 
reproduces the models in the theoretical HR diagram for comparison. 


gram above the standard deuterium birth line therefore 
represents a “forbidden region”, in which stars with low 
spot coverage may only reside for a short time. By con¬ 
trast, heavily spotted stars naturally occur in this loca¬ 
tion due to their intrinsically lower average surface tem¬ 
peratures. For example, the spotted 3 Myr isochrone in 
Fig. piies entirely above the un-spotted birth line below 
^ O.5M0. Such stars persist in the forbidden region for 
several tens of Myr, given their slow contraction rates, 
and should therefore dominate the population coolward 
of the standard deuterium birth line. This feature pro¬ 
vides a direct observational probe of the prevalence of 
spotted stars in young clusters. 


4.2. Spot effects in the color-magnitude diagram 


In 14.1 we drew T^e and luminosity directly from our 
models to predict HR diagram shifts. While valid if 
and L are known, this approach neglects an important 
fact: spotted stars are inhomogeneously hot across their 
surface. This will have implications for the spectrum of 
emission from the photosphere, and may introduce color 
anomalies relative to pristine stars of equal Teff. To in¬ 
vestigate this issue, we transform our models into color 
space using a two-zone model, where the emission is the 
area-weighted composite of the SEDs from the surfaces 
of temperature Tamb and Tgpot, whose values are derived 
from the model Teff using Eq. This transformatio n is 
accomplishe d by using the bolometric corrections of An 
et al. (20071. Colors are rarely used to interpret pre-Ms 


stars as numerous complicating factors (circum-stellar 
discs, heavy extinction, etc.) render them un-reliable. 
We therefore advance our models to 125 Myr for compar¬ 
ison with data in the Pleiades, but note that analogous 
effects occur on the pre-MS. 

The top left of Eig. shows the results for My vs. 
V — K. The red line represents un-spotted models from 
O.5-1.OM0, and the dark blue line shows our spotted 
models, transformed using the two-zone approach (xgpot 
= 0.8, /spot = 0.5). The spotted models appear redward 
of the standard models, though the shift is not as extreme 
as in the theoretical HR diagram (shown in the bottom 
right). This is because the un-spotted surface temper¬ 
ature is higher than the model Tgff, thus enhancing the 
short wavelength emission and producing bluer colors. 


These models ar e compared to Pleiades data from Ka- 
([20M ) sup plemented by rotation data from 


mai et al. 


Hartman et al.| (|2010|). Pleiads rotating faster than a 
two day period are snown as open circles, and slower 
stars are shown as black dots. The slow rotators gen¬ 
erally agree with our un-spotted models, but the fastest 
rotators are typically more red by ^0.1-0.2 magnitudes. 
As our spotted models also appear more red in this color 
band, this supports the interpretation that spots are af- 


Stauffer et al. 

2003 

Jackson & Jeffries 

2014a 


we u nder-predict the magnitud e of the shift. 'The mod¬ 
els of Jackson & Jeffries (2014a I found that spots induce 
a blue-ward shift in this plane, in contrast to the red- 
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Figure 5. The effect of spot filling factor and temperature on the color-magnitude locations of spotted models with the same spot blocking 
parameter aspot- Darker spots cause models to appear bluer, as the hot, un-spotted photosphere progressively dominates the total emission. 
The bottom right panel shows surfaces of constant flux redistribution in the plane of fspot vs. Xspot, and the colored dots represent the 
different models shown in the other three panels. 


ward shift we find. This discrepancy may be related to 
differing choices of the spot temperature contrast (see 
below). 

Next, we present isochrones in two additional color co¬ 
ordinates: My vs. V — Ic, and My vs. B — V. As can be 
seen, the spotted models appear bluer in shorter wave¬ 
length colors. The spotted isochrone lies nearly on top of 
the un-spotted isochrone in V — Ic, and completely blue- 
ward of it in B — V, similar to the behavior of the rapidly 
rotating data. This is due to the progressively stronger 
sampling of the Wein tail as one considers shorter wave¬ 
lengths. It has long been known that the K dwarfs of 
the Pleiades are blue in B — V, and red in V — K, com¬ 


pared to i sochrones calibrated on older stars. [Stauffer et 
al. (20031 attributed this behavior to spots and plages 


Brought on by magnetic activity, which respectively en¬ 
hance the long and short wavelength emission. Our mod¬ 
els suggest that qualitatively similar behavior is possible 
without plages, as the extra blue light can result from 
a slight increase in the temperature of the un-spotted 
photosphere. In reality, there are several effects we have 
neglected in our models which impact the detailed lo¬ 
cation of CMD isochrones, such as plages, faculae, and 
other activity-related phenomena. A more complete and 
detailed model is therefore necessary to make definitive 
quantitative statements. 

These two-zone models were calculated using the s ame 
combination of fspot and Xgpot- As discussed in ([2.2 


these parameters are degenerate from a modeling per¬ 


spective, since cooler spots with a lower surface area can 
redistribute the same amount of flux as warmer spots 
with a higher surface area. This degeneracy is CTaph- 
ically represented in the bottom right of Fig. ^ the 
curved lines denote surfaces of constant flux redistribu¬ 
tion (our base case implies aspot ~ 0.7). However this 
modeling degeneracy does not require that the param¬ 
eters be observationally degenerate, as the spectrum of 
photospheric emission depends on the temperatures and 
relative fractional coverage of the two zones. 

We therefore present CMD isochrones calculated with 
different degenerate combinations of Xgpot and fspot in 
Fig. 0 Each panel shows our un-spotted isochrones 
in red, and isochrones derived through the two-zone ap¬ 
proach for four degenerate choices of Xspot and fspot (see 
the colored dots in the bottom right panel). In each case, 
we see that lower spot temperatures (lower Xgpot) give 
rise to bluer emission, since the decreasing flux of the 
spots leads to increasing dominance of the global emis¬ 
sion by the hot ambient surface. This effect saturates 
at sufficiently low Xspot and asymptotes to totally black 
spots. The saturation value is a strong function of the 
color bands: for example, Xspot = 0.6 models are nearly 
indistinguishable from black spots in B — V, while they 
remain separated in V — Ic and V — K. In the limit¬ 
ing case of black spots, the entire flux budget originates 
from the ambient regions, leading to bluer emission in all 
bands. This behavior is consistent with the findings of 
Spruit & Weiss (19861, who considered the effect on col- 
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Figure 6. The impact of spots on the accuracy of cluster ages obtained through the lithium depletion boundary technique. Left: Lines 
of constant spot filling factor denote the age at which 95% of the primordial Li is destroyed from models with a given luminosity. Spotted 
models lie above the un-spotted case, demonstrating that spots delay the onset of Li destruction, thus biasing determinations towards 
younger ages. Right: The fractional error on the age determination as a function of true age, for three different filling factors. The errors 
are approximately half of the spot blocking parameter aspot (0.3, 0.2, and 0.1 respectively for purple, blue, and red) up to ~60 Myr, and 
drop off thereafter. 


ors of black spots (see their fig. 2), but did not explore 
the influence of temperature contrast. 

The fact that the observed emission of spotted stars is 
a composite of two zones may hold crucial implications 
for the spectra of young stars. Ionization fractions will 
be non-uniform across the surface, the strength of Zee- 
man broadening and splitting will differ wildly from loca¬ 
tion to location, and the cool spotted regions may permit 
molecular species which are too fragile for the hot photo¬ 
sphere. Given these effects, it is clear that the observed 
spectrum will be different than a pristine star of equiv¬ 
alent Teff, which may complicate the inference of spec¬ 
troscopic stellar parameters and metal content. Further¬ 
more, these isochrones show that multi-band photometry 
holds great value as a tool for breaking the degeneracy 
between /spot and Xgpot, which limits the conclusions of 
many observational studies of starspots. For example, a 
given spot modulation signal may increase due to greater 
(non-axisymmetric) spot coverage, or due to cooler spots. 
Combining spot modulation with isochrone comparisons 
may help to constrain the relative contribution of spot 
area and spot temperature. While we have exhibited the 
importance of this effect, a full exploration of requires 
additional important stellar physics as described above, 
and is therefore outside our scope. 

5. EFFECT ON LITHIUM DEPLETION AND CLUSTER 
DATING 

Lithium has long been used as an indicator of youth in 
late-type stars. Li is destroyed at T^i ^ 2.5 x 10® K for 
typical stellar densities, and consequently depletes from 
the surface of stars when the base of the convective en¬ 
velope reaches this temperature on the pre-MS. For stars 
more massive than ^O.SMq, the convection zone retreats 
and cools below T^i near the end of the pre-MS, thus 
freezing out the surface Li abundance and imprinting a 
characteristic mass-dependent pattern at the ZAMS. For 
stars between O.SMq and the sub-stellar boundary, Li 
in the outer layers will be completely destroyed before 
reaching the MS. Li is also destroyed for some objects 


below the sub-stellar boundary (O.O8M0< M <O.O6M0), 
but lower mass brown dwarfs never reach T^i in their in¬ 
terior, and thus do not destroy their Li. As the rate of 
burning is extremely sensitive to temperature (oc 
any process that impacts the internal structure of stars 
may have implications for Li abundance predictions. In 
this section, we explore the impact of spots on Li destruc¬ 
tion in two mass regimes: fully convective stars, which 
have been widely employed to measure cluster ages, and 
FGK dwarfs, whose depletion properties have long been 
at odds with standard theoretical predictions. 

Perhaps the most important utility of Li is the mea¬ 
surement of young cluster ages throu gh the lithium de¬ 


pletion boundary (L DB) technique (Basri et al. 1996 


Bildsten et al. 1997). When the central temperature o: 


a fully convective star reaches on the pre-MS, Li is 
destroyed throughout the star in just a few Myr. The 
rate of Hayashi contraction, and thus the precise age 
at which Li destruction occurs, depends strongly on the 
masses of individual stars. As a consequence, there ex¬ 
ists at a given age a critical mass below which stars re¬ 
tain their primordial abundance, and above which Li has 
been annihilated. Identification of this transition mass, 
and the subsequent inference of cluster age by compar¬ 
ison with stellar models, comprises the LDB technique. 
This method has been favored by many authors, due to 
its relative insensitivity to errors in the physical input s 
of stellar evolution codes (~ 3 — 8% |Burke et al.||2004|). 
However, Jackson & Jeffries (2014b) recently showed that 
because spots reduce the rate of pre-MS contraction (sim¬ 
ilar to our findings), the age at which T^i is reached in 
the interior of fully convective stars increases, introduc¬ 
ing a bias into LDB ages inferred from spotted stars. 
In the following, we perform a similar anal ysis w ith our 
evolution code, and compare the results in §6.2| 

To quantify this effect, we present models of M < 
0.5Mq in the left panel of Fig. The lines in this 
figure represents the age at which 95% of the primordial 
Li abundance is destroyed as a function of bolometric lu- 
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Figure 7. The impact of pre-MS spots on Li destruction in par¬ 
tially convective stars. Higher coverage results in a lower pre-MS 
depletion rate, and thus a higher abundance at the ZAMS. The 
resulting range and temperature dependence is similar to the Li 
pattern observed in the Pleiades (grey circles). 


minosity in un-spotted (black line/squares) and spotted 
(magenta lines/squares) models. The fractional delay in 
the onset of Li destruction is greater than the fractional 
luminosity reduction, implying that derived ages will be 
under-estimates. The right panel of Fig. shows the 
fractional age error incurred by using un-spotted evolu¬ 
tionary models to measure the lithium depletion bound¬ 
ary from stars with three different spot filling factors, 
as a function of the true age. The measurement pre¬ 
cision remains roughly constant from ~15 to 60 Myr, 
and decreases sharply thereafter, due to the pro gressively 
stronger L reduction in lower mass stars (j ]3.1[ ). For our 
maximal spot coverage, the fractional age error begins 
around 13% at 15 Myr, peaks at ~ 15% at 30 Myr, and 
decreases below 5% after 100 Myr. 

These numbers indicate errors for a uniformly spotted 
population. If spot properties vary from star to star, 
the effect will be to eliminate the one-to-one mapping 
between mass and Li destruction age, thus blurring out 
the otherwise sharp depletion boundary. The strength 


of activity may also depend on sp ectral type (Mohanty 
fc Basrij2003[ [Schmidt et al.|2015 ), particularly near the 
sub-stellar boundary (Sp'L ^ M8) in which case different 
values of /spot should be considered for different masses of 
the depletion boundary. LDB ages have long been in con¬ 
flict with ages derived from isochrones, in the sense that 
LDB ages are systematically higher than MS isochrone 
fitting, and the spot-induced errors serve to exacerbate 
this conflict. We note however that the two techniques 
can be brought into better agreement by the inclusion of 


convective overshooting on the upper MS (e.g. Soderblom 
eFaI][M4l ). 


IN ext, we evaluate the impact of spots on Li destruction 
in higher mass stars. As discussed above, stars which de¬ 
velop radiative cores cease surface Li depletion when the 
convective envelope retreats and cools near the end of 
the pre-MS. This occurs at earlier ages for higher mass 


stars, so the amount of Li preserved at the ZAMS in¬ 
creases with effective temperature. Standard theoretical 
models anticipate a strong mass tren d with little scat¬ 
ter at the ZAMS ( Pinsonneault||1997 ). Contrary to this 
picture, large (> lUx j Li dispersions at fixed Teg have 
been identified in the K-dwarfs of several young cluster s 
(e.g. Soderblom et al. 1993 Balachandran et al. 2011). 
Such dispersions further correlate with rotation, so that 
the fastes t spinn ing stars tend to be the most Li-rich. We 
argued in |SP14| and |SP15| that if some rotation-correlated 
mechanism inflated stars relative to standard predictions 
on the pre-MS, then a rotation-correlated Li dispersion 
would naturally emerge in open clusters around 10-15 
Myr. A dispersion in Li at the ZAMS thus represents 
a fossil record of non-uniform stellar parameters among 
pre-MS stars during the epoch of Li burning. We have 
now implemented a specific m echa nism which predicts 
radius inflation on the pre-MS ((3.1), and which is known 
to correlate with rotation in some stars, so we examine 
its effect on Li depletion in K-dwarfs. 

Fig. Ill shows the Li abundances of O.6-1.2M0 stellar 
models at 125 Myr (the age of the Pleiades). Each line 
represents a different fspot, given by the key in Fig. P 
{Xspot = 0.8 in each), and each dotted line connects mom 
els of equal mass. As anticipated, different filling factors 
lead to different Li abundances at the ZAMS, with the 
dispersion surpassing a factor of ten at some masses. We 
have also plotted A(Li), and rotation infor mation 
for a sample of stars in the Pleiades , ass embled in|SP15 
from d ata in Soderblom et al. (|1993 ) and Hartman et al. 
(2010). Point size indicates rotation period, from which 
the Li-rotation correlation can be readily identified be¬ 
low ~ 5000 K. As can be seen, the abundance dispersion 
in these data roughly corresponds to models in the range 
fspot= 0.0-0.45 {aspot = 1.0 — 0.73). This analysis sug¬ 
gests that if rotation correlates with spottiness on the 
pre-MS as it does on the MS, spot-induced radius in¬ 
flation naturally produ ces the Li-rotation correlation in 
this cluster (see |SP15| for a lengthy discussion). This 
in turn implies an early (^3 — 15 Myr) origin for the 
Pleiades Li dispersion, a hypothesis which may be tested 
by obtaining Li abundances for clusters of age ~10-15 
Myr. Finally, while many stars are more depleted than 
the un-spotted models above 5000 K, th is disc repancy 
could be explained by rotational mixing (SP15I, errors 
in the input physics of the models (SP14), or could indi¬ 
cate that the Pleiades was born with a sub-solar initial 
Li abundance (A(Li) < 3.31). 

6. DISCUSSION 


We have implemented a simple spot model into our 
stellar evolution code, and assessed its effect on stars dur¬ 
ing the pre-MS and on the ZAMS. We have particularly 
focused on the impact of spots on, 1) stellar structure 
as a function of mass; 2) the isochrone-derived masses 
and ages of pre-MS stars; 3) The colors of spotted stars 
near the ZAMS; 4) the depletion of lithium from stel¬ 
lar envelopes. In this section, we synthesize our results 
and consider some implications, compare our results with 
other works who consider magnetic models, and discuss 
some outstanding issues in both the physical interpreta¬ 
tion and modeling of spotted stars. 


6.1. A radius dispersion during the pre-main sequence 
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We have shown that if a range of spot properties (i.e. 
surface magnetic fields) exist at fixed mass and age dur¬ 
ing the pre-main sequence, the result is non-uniform stel¬ 
lar parameters among the members of young open clus¬ 
ters. But does this scenario accord with observational 
evidence? 

In the case of lithium, we find a conspicuous spread 
in abundances at fixed temperature among the G- and 
K-dwarfs of open clusters shortly after the age of pre-MS 
Li burning. This sp read cannot be ex plained by atmo¬ 


spheric effects (e.g. King et al. 2010), by Li depletion 
in non-magnetic models ([SBldl), or by the inclusion of 


rotational mixing (SP15[), but naturally emerges from a 


range of spot properties (i.e. a range of radii) at fixed 
mass and age during the pre-MS. Furthermore, our pre¬ 
dictions are consistent both with the radius data for ac¬ 
tive eclipsing binaries, and with the anomalous colors of 
young active stars in the Pleiades, suggesting that our 
models can reproduce observational results with realistic 
spot properties, and that magnetic activity continues to 
inflate rapidly rotating stars on the ZAMS and beyond. 
Collectively, these results indicate that the differential 
magnetic properties of stars (correlated with rotation) 
lead to a dispersion in radius at hxed mass and age, 
which emerges on the early pre-MS and persists on the 
MS while stars remain rapidly rotating. 


radius dispersion in pre-MS 1 

Cottaar et al. 

2014), and 

young MS ( 

Jackson et al.|2009 

Jackson &: Jet 

irles|2U14al 


appear to be m conhict with 
the uniform saturation of magnetic proxies among young 
stars, which could be interpreted to mean that spot prop¬ 
erties are identical at fixed mass on the pre-MS. However, 
we note that saturation is merely an empirical statement 
about chromospheric and coronal heating, and that a di¬ 
rect correlat ion w ith spot properties has not yet been 


proven (see (6.3.2). Alternatively, saturation may indi- 
cate that oth^pEeno menon, such as mass accretion (e.g. 


Littlefair et al. 2011), play a role in generating the pre- 
MS radius disper sion n eeded to explain the Li data (see 


the discussion in SP15), though this scenario still relies 
on activity to explain the radii of eclipsing binaries and 
the colors of rapidly rotating Pleiads. 

A radius dispersion in pre-MS clusters have several sig¬ 
nificant implications. First, we have shown that spots 
make stars appear cooler, and as a result heavily spot¬ 
ted stars coincide on the HR-diagram with younger , les s 


massive stars which are weakly spotted or pristine ((4.1). 


The range in masses at fixed HR diagram location can 
reach a factor of two, while the range in ages can be sig¬ 
nificantly higher (up to ^ 3x above O.SM©, and larger 
for lower masses). Consequently, pre-MS cluster ages 
inferred by matching isochrones to the locus of HR di¬ 
agram data will be erroneously low. The magnitude of 
the age error may match the above-quoted correction fac¬ 
tors if stars are uniformly spotted, or may be somewhat 
less if a range of spot hlling factors blurs out an other¬ 
wise sharper stellar locus. Furthermore, if stars of equal 
age have differential spot properties, the inferred ages 
can differ substantially, leading to the appearance of an 
intra-cluster age spread. Appreciating the magnitude of 
this effect is important for inferring the true age spreads 
within young clusters. 

As noted in the lithium depletion boundary tech¬ 


nique also systematically under-estimates the ages of 
clusters whose members are heavily spotted, by up to 
15% in the most extreme case we considered. As many 
previous young cluster studies do not consider these mag¬ 
netic effects, our findings necessitate an upward revision 
of the age scale for clusters between 1 and 100 Myr. The 
magnitude of this revision is difficult to anticipate, as 
the required corrections for each individual study depend 
strongly on the true age of the cluster and on the mass 
range probed, as well as the distribution of spot prop¬ 
erties among stars in each cluster. A full re-analysis of 
young cluster CMDs in the context of magnetic models 
is beyond our current purpose, but would undoubtedly 
have important implications for, among other things, the 
duration of circum-stellar disc lifetimes and the accre¬ 
tion timescales of proto-stars, both of which are inferred 
through cluster ages. 

Second, Fig. shows that among spotted stars, the 
least massive require the largest correction factors. Con¬ 
sequently, in a mono-aged cluster the lowest mass objects 
will appear to have the youngest isochronal ages, and the 
hottest stars will appear to have the oldest isochronal 
ages. Thus in the presence of spots, clusters should show 
a trend in inferred age with temperature. This precise 
behavior was note d in the young /3 Pic association by 


Malo et al. (2014): stars with Tes> 3500 K appeared 
to be 15 Myr on average, while cooler stars appear to 
be 4.5 Myr on average, both of which are yo unger than 
the generall y accepted age of ^ 21 ± 4 Myr (Soderblom 


et al. 2014). They showed that this discrepancy could 
be arneliorated through the inclusion of strong magnetic 
fields in the cooler stars, and weaker fields in the hot¬ 
ter stars. Alternatively, consistent spot properties at all 
masses will naturally impact the age determinations of 
lower mass stars to a greater degree, and thus can pro¬ 
duce just such a disagreement. The relative isochronal 
ages of cool and warm stars depends on the relative spot 
properties of the two populations, and therefore provides 
an empirical constraint on the surface magnetic fields of 
pre-MS stars. 

Finally, the systematic under-estimation of the masses 
of stars will bias measurements of the initial mass func¬ 
tion (IMF) which employ pre-MS clusters. In particular, 
the average low mass star will require a significant up¬ 
ward revision of its mass, ultimately implying a more top 
heavy IMF, and a lower number of brown dwa rfs. The 
magnitude o f this effect has been discussed in |Stassun| 


et al.| (2014), and we refer the reader there for a more 
detailed account. 

6.2. Literature comparison 

Heavy spot coverage can enhance stellar radii by up 
to ~10% during the pre-MS at fixed age, due in part to 
the slower contraction and in part to the impact of spots 
on the pressure conditions at the photosphere. Once on 
the MS, the peak radius anomalies are mass dependent, 
with partially-convective stars generally more inflated 
than fully-convective stars. Thi s pattern is unsurpris - 
ingly similar to the findings from Spruit & Weiss (1986), 
with the exception that our models predict a decrease m 
AR at hxed sp ot properties between O.4-O.8M0, while 
[Spruit &: Weiss| predict flat Ai? in this range. The origin 
of this discrepancy is unclear, but is probably related to 
different input physics which cause a substantial shallow- 
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ing of the mass-radius relationship around O.TMq. 

Assuming equipartition between the magnetic and 
gas pressures in the spotted regions (i.e. Pspot.mag = 
Pspot,gas) and pressure equilibrium between the spots 
and the photosphere, our models imply magnetic field 
strengths of 2.2, 1.6, and 1.1 kG at the ZAMS for 0.2, 
0.5, and I.OMq, respectively. This is generally consistent 
with the 1- 3 kG range reported fo r empirical measure¬ 
ments (e.g. Cranmer & Saar 2011). Gombined with our 
spot coverage factor of 0.5, the surface averaged fields 
<!?«/» > = 0.6-1.1 kG. This value i s similar to the find¬ 
ings of Mullan & MacDonald (2001|, who require surface 
strengths of less than 1 kG to inflate fully convective stars 
by ~3-4%. They are however some what lower than the 
2.0-6.0 kG s urface fields required by [Feiden fc Ghaboyer] 
(2013 2014) to match the observed radii, T^gs, and lu¬ 
minosities of several eclipsing binaries. This quantita¬ 
tive discrepancy is not surprising, as their studies make 
significantly different assumptions about magnetic field 
morphology. We note however that the mass dependence 
they find is qualitatively similar to ours; stronger mag¬ 
netic fields are required to inflate fully convective stars 
relative to solar type stars, consistent both with the pat¬ 
tern revealed in Fig. and the estimates given 

above. This suggests that greater inflation at higher mass 
is a generic prediction of magnetic stellar models. 

In ^ we found that spot-induced errors on cluster 
ages derived from the LDB technique are large com¬ 
pared to systematic uncert ainties in the mode l inputs 
(^15%, compared to <8%; |Burke et al.||2004D. How¬ 
ever, the errors we derived are somewhat lower than the 


^20% age errors reported by Jackson & Jeffries (2014b), 
who considered the impact of similar spot coverage on 
lithium depletion in their polytropic models. These re¬ 
sults are in reasonably good agreement, considering the 
numerous physical inputs required in stellar modeling, 
but it is interesting to consider the source of this dis¬ 
crepancy. In a polytrope model, the internal gradient 
V is assumed to be Vad throughout the entirety of con¬ 
vective regions. While this is a good assumption for the 
vast majority of the interior of fully-convective stars, it 
is less exact in the thin layer at the surface where the 


gradient beco mes super-adiabatic (e.g. Kippenhahn & 
Weigert 1990). As it happens, starspots chiefly reside 


near this super-adiabatic layer, and so their effect on the 
structure of the star will be dependent on the adopted 
physics of this region. As discussed in ^ changing the 
temperature at the boundary of a convective zone leads 
to a linear shift in P throughout the region. Starspots 
are placed at this o uter bo undary in the polytropes of 
Jackson & Jeffries (2014bI, so the resulting L shift in 
the center is approximately equal to the L blocked at 
the surface. In our standard stellar models, the temper¬ 
ature shift first propagates through the super-adiabatic 
layer before impacting the convective envelope, leading 
to a somewhat higher central value of L. This in turn 
impacts the Li depletion age as a function of mass, and 
thus the inferred ages of clusters. However, we note that 
the different quantitative results may derive simply f rom 
different choices of input physics (Burke et al.||2004' I. 


6.3. Observational and physical concerns 

Our approach in this paper has been to adopt a fixed 
set of spot parameters (i.e. fspot = 0.5, Xspot = 0.8 re¬ 



Figure 8. The impact of spot depth on radius inflation. The black 
and red dashed lines represent shallow and deep spot treatments, 
respectively. The top shows absolute anomalies, and the bottom 
shows the difference between the anomalies as a function of mass. 
The difference reaches at most ~1.5%, and thus impacts detailed 
calculations, but not the conclusions of this paper. 

gardless of age), and derive their effect on stellar prop¬ 
erties. While such limits are useful, the ideal analysis 
would tie the temperature, covering fraction, and distri¬ 
bution of spots to the (evolving) structure and properties 
of their host stars. This physical approach would permit 
a realistic population model of spotted stars in young 
clusters, which is desirable for understanding the effect 
of spots on observables during the pre-MS and MS. Given 
the present state of the field, such a model awaits supe¬ 
rior understanding of the underlying mechanisms govern¬ 
ing the behavior of stellar magnetism. We highlight in 
this section some key outstanding questions whose reso¬ 
lution will improve the capacity of stellar theory to model 
magnetic fields in young clusters. 


6.3.1. Spot depth 

Spots extend to a finite depth within the envelopes 
of stars, but exactly how deep they penetrate, and the 
details of their subsurface morphology, remain a topic 
of debate (e.g. Rempel & Schlichenmaierl 20111. In our 
models, we have assumed that 'lamb — d'spot is constant 
as a function of depth, a criterion which leads to shallow 
spot^ To assess how this assumption impacts our re¬ 
sults, we compare our base case to models which assume 
a constant temperature contrast with depth {Tamb/Tspot 
= constant). This criterion implies very deep spots, since 
the temperature contrast at the surface is maintained at 
the base of the surface convection zone. 

The comparison appears in Fig. The top panel 
shows the radius anomalies resulting from the two depth 
treatments at the ZAMS, and the bottom panel reports 
the difference between the curves. Deep spots lead to 
greater radius inflation with a maximum difference of 
^1.4%. This is a small, but non-negligible effect relative 


® As an example, a surface spot contrast of Xspot = 0.8 in a IMq 
model reduces to less than a 2% effect at about 3 pressure scale 
heights below the photosphere, being in the outer 10~® and 10“® 
in mass and radius, respectively. 
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to the overall radius anomaly (~10%). This suggests that 
the global impact of spots is chiefly due to their effect 
on the very outer regions of the star, but is also moder¬ 
ately sensitive to flux redistribution in the envelope. Our 
constant spot contrast is an extreme limiting case, as it 
would imply an extremely strong magnetic field (tens to 
hundreds of MG), and a strong magnetic field gradient, 
near the base of the convection zone. We therefore ex¬ 
pect the true difference to be significantly less than seen 
in Fig. 1^ Despite the marginal differences, shallow and 
deep spots produce quantitatively similar results, and we 
therefore conclude that depth effects do not impact our 
results. 


6.3.2. Magnetic saturation 

A key outstanding question in the study of stellar mag¬ 
netism is the physics underlying magnetic saturation. 
Numerous studies have determined that once stars are 
rotating faster than a critical Rossby num ber (^^0.1; e.g. 
Noyes et al. 1984 Pizzolato et al. 2003), classical mag- 
netic field proxies cease to increase with faster rotation. 
This observation suggests a saturation in the strength 
of magnetic torques and chromospheric/coronal heating, 
which are the physical processes that generate the ob¬ 
served proxies. The cause of this saturation remains un¬ 
clear, but is probably related to a change in the proper¬ 
ties of the stellar dynamo when the rotation period be¬ 
comes substantially shorter than the convective overturn 
timescale (i.e. low Rossby number). 

Observations of young clusters indicate that the vast 
majority of pre-MS stars are in the satu rated regime 
(Preibisch et al. 2005 Argiroffi et al. 2014). If all mag- 
netic effects, such as spots, saturate along with the 
classical proxies, then one would expect uniform spot- 
properties at fixed mass on the pre-MS, and thus no 
dispersion in stellar properties due to magnetic activity. 
However, there is at present no convincing demonstra¬ 
tion that spots saturate along with magnetic proxies, 
and some studies have in fact found evidence that the 
amplitude of spot mod ulation varies substantially in the 
saturate d regime (e.g. O’dell et al. 1995 McQuillan et 
al.||2014|), perhaps suggesting disparate spot properties, 
furthermore, a uniform spot coverage does not accord 
with the evidence from Li observations, which indicate a 
rotation-correlated dispersion in radii at fixed mass and 
age during the pre-MS (l|^. Uncovering the spot prop¬ 
erties of pre-MS stars has been hampered by observa¬ 
tional issues, primarily the substantial degeneracies be¬ 
tween filling factor, the surface distribution of s pots, and 
their temperatu re relative to the photosphere (Jackson 
& Jeffries 2012). Nevertheless, the connection (or lack 
thereof) between saturated proxies and spot properties 
is an essential component of a complete theory of stellar 
magnetic activity, and must be understood in order to 
appreciate the impact of spots on the pre-MS. 


6.3.3. Plages 

Another confounding effect is the presence of plage re¬ 
gions in stellar photospheres. Plages are the hot analogs 
of starspots, and radiate light at higher frequencies than 
the un-spotted surface. This enhanced emission may 
counteract the total flux suppression caused by spots, 
and thus reduce the magnitude of the perturbations to 
stellar parameters that we have outlined in l|3][^ The 


size of this correction relies on the difference between the 
plage enhancement and spot decrement to the total irra- 
diance, which is likely dependent on mass and rotation 
among other parameters. Furthermore, plage regions 
provide a source of anomalous blue emission, which has 


significant effec ts for the SEDs of active stars (e.g. Stauf- 
fer et al.||2003|. These regions are poorly understood at 


present, so a deeper understanding of their origin, and 
their correlations with spots on pre-MS and ZAMS stars, 
is necessary. 


6.3.4. Spot evolution 

Finally, in our models we have assumed that spots ap¬ 
pear before the deuterium birth line, and persist at the 
same magnitude for the entirety of the approach to the 
MS. In reality, spot properties likely change over contrac¬ 
tion timescales as the rotation rate and convective prop¬ 
erties of stars evolve, as well as on shorter timescales due 
to the appearance and dissolution of individual spots. A 
fixed filling factor can approximate the time-averaged ef¬ 
fect of spots, but may miss crucial physics occurring on 
sub-thermal timescales. Furthermore, we have assumed a 
single temperature contrast between the spotted and am¬ 
bient regions. In fact , spot temperatur es appear to cor¬ 
relate with mass (e.g. Berdyugina|2005 ), and may evolve 
during the pre-MS due to the changing pressure, density, 
and rotation rate of contracting stars. A superior treat¬ 
ment of these dynamic quantities is necessary for robust 
spot predictions. 

The heliographical distribution of spots may also play 
a role in determining the efficiency of energy redistribu¬ 
tion. For instance, if a large number of small spots are 
evenly distributed throughout the photosphere, then the 
entire un-spotted surface will be in close proximity to 
magnetic regions. By contrast, if the star is populated 
by a few large spots, then certain ambient regions will 
be much farther from the closest magnetic region. The 
average distance of the photosphere from spots could 
impact the efficiency and uniformity of flux redistribu¬ 
tion, and might therefore change the structural effect of 
starspots. More sophisticated (2-D or 3-D) spotted mod¬ 
els may help shed light on the role of spot distribution 
on surface and sub-surface structure. 


7. SUMMARY AND CONCLUSIONS 

Through the inclusion of spots in our stellar models, 
we have explored the impact of a non-uniform surface 
on the radii, Tgfj, luminosity, inferred masses and ages, 
colors, and Li abundances of pre-MS and ZAMS stars. 

We find that a 50% coverage of 80% temperature con¬ 
trast spots (an effective blocking area of 30% black spots) 
increases the radii of all stars on the pre-MS by 8-12%, 
and by a mass-dependent amount on the MS, which in¬ 
creases from ^3% in fully-convective stars to ~10% at 
1.2M0. This behavior is a product of both the impact 
of spots on the central luminosity of stars, and on the 
fine structure of the mass-radius relationship. Moreover, 
these anomalies are similar in magnitude to claimed ra¬ 
dius anomalies in detached eclipsing binary systems, sug¬ 
gesting that spots may play a role in the well known 
radius inflation of young, active low mass stars. 

In our models, these radius anomalies are accompanied 
by corresponding decrements to the luminosity, which 
are large in fully-convective stars and small in stars with 
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thin envelopes, and to the Teff, which are weakly sensitive 
to mass. These perturbations are physical, and displace 
stars from their fiducial location in the HR diagram. This 
effect has several implications: 

1) Masses derived from spotted stars using pre-MS 
isochrones will be systematically low by up to 2x at 
3 Myr, and ages will be systematically young by up to 
10X for the lowest mass stars. We quantify the result¬ 
ing errors as a function of mass and age, and present 
correction factors for spotted stars. 

2) Spotted stars will be preferentially displaced to¬ 
wards the upper right of the standard HR diagram; in 
particular, they may naturally appear above the un¬ 
spotted deuterium birth line. Detecting spot modula¬ 
tion in this region of the HR diagram therefore presents 
a direct test of the prevalence of spots on pre-MS stars. 

3) A range of masses and ages may be present at fixed 
HR-diagram location, and low mass stars will generally 
appear younger at fixed age compared to higher mass 
stars. Interpreters should be cautious about relaying the 
fundamental parameters of stars without considering the 
impact of activity. 

4) The observed colors of stars will reflect an area- 
weighted sum of the emission from the hot and cold re¬ 
gions. We find that spotted stars appear bluer than un¬ 
spotted stars in short-wavelength bandpasses, and red¬ 
der in long-wavelength bandpasses, consistent with ob¬ 
servations of (spotted) stars in the Pleiades and other 
young clusters. Furthermore, the extent of the blue-ward 
shift depends strongly on the average temperature of the 
spots, which may help to break observational degenera¬ 
cies between total spot coverage, surface distribution, 
and temperature. 

Finally, we examined the impact of spots on pre-MS 
lithium depletion in two regimes: fully-convective stars, 
and in G and K dwarfs. For the former category, the 
main interest is measuring the ages of young clusters 
through the lithium depletion boundary technique. We 
found that an effective blocking area of 30% can in¬ 
crease the inferred age by ~ 10-15% relative to the true 
age. This bias exacerbates the tension between LDB 
ages and MS turn-off ages in young clusters such as the 
Pleiades. The errors we derive are less severe than the 


results of Jackson & Jeffries (2014b), with the difference 
potentially stemming from the treatment of the super- 
adiabatic region in the two modeling approaches. For G 
and K dwarfs, we find that spots inhibit the destruction 
of Li during the pre-MS by altering the temperature at 
the base of the convection zone. A range of spot filling 
factors in our pre-MS models leads to varying degrees of 
Li destruction and thus a dispersion at fixed mass and 
temperature on the ZAMS, consistent with observational 
studies of young clusters. This implies that clusters dur¬ 
ing the epoch of Li destruction (~ 5 — 15 Myr) are dif¬ 
ferentially affected by magnetism, imprinting a radius 
dispersion of ~ 5 — 10% at fixed mass. 
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